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Predicted gamma-ray image of SN 1006 due to inverse 
Compton emission 
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ABSTRACT 

We propose a method to synthesize the inverse Compton (IC) gamma-ray image of 
a supernova remnant starting from the radio (or hard X-ray) map and using results 
of the spatially resolved X-ray spectral analysis. The method is successfully applied 
to SN 1006. We found that synthesized IC gamma-ray images of SN 1006 show mor- 
phology in nice agreement with that reported by the H.E.S.S. collaboration. The good 
correlation found between the observed very-high energy gamma-ray and X-ray /radio 
appearance can be considered as an evidence that the gamma-ray emission of SN 1006 
observed by H.E.S.S. is leptonic in origin, though the hadronic origin may not be 
excluded. 

Key words: ISM: supernova remnants - shock waves - ISM: cosmic rays - radiation 
mechanisms: non-thermal - acceleration of particles - ISM: individual objects: SN 1006 
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1 INTRODUCTION 

H.E.S.S. and MAGIC projects open the era of the very-high 
energy (VHE) 7-ray astronomy in the sense that they can 
provide for the first time detailed images of various astro- 
physical objects. VHE 7-ray emission from supernova rem- 
nants (SNRs) is one of the key components in investiga- 
tion of the processes around strong nonrelativistic shocks, 
namely, dynamics and structure of the shock itself, mag- 
netic field behavior, and microphysics of charged particles 
including their injection and acceleration. Observations in 
this energy domain demonstrate that charged particles are 
really accelerated in SNRs up to the highest energies ob- 
served in galactic cosmic rays. 

SNRs with kno wn TeV 7-ray images includ e 
RX J171 3 7-3946 (lAharonian et all 1200a l2007aT ). 
Vela Jr. (lAharonian et all l2005cL l20Q7bh. SN 1006 
(iNaumann-Godo et al.1 12009). R CW86 (lAharonian et~al 
091), IC443 ([Albert et al.1 l2007h. W28 (lAharonian et~ 
08ah. CTB 37B (lAharonian et all l2008bh, GO.9+0.1 
(jAharonian et al.1 l2005aT ). MSH 15-52 (|Aharonian et al.1 



2005b). Pulsar-wind nebulae could be the origin of 7-rays 
in GO.9+0.1 and MSH 15-52 while the particles accelerated 
at the forward shock of SNRs are likely to be responsible 
for emission in the other shell-type SNRs. 

The broad-band analysis of the spectra from SNRs 
is useful to set constraints on model parameters but it 
still leave open the nature of VHE 7-ray flux either as 
leptonic or as hadronic in origin (e.g. RX J1713. 7-3946: 
lAharonian et al.l l2006: Berezhko & Volk 2006). The analysis 
of the spatial distribution of 7-ray emission is an additional 
important channel of the experimental information. 

Hadronic 7-rays arise at the location of the target pro- 
tons. Rather large density of target protons - as e.g. in 
molecular clouds - is the condition for the effective hadronic 
emission in SNRs with high TeV 7-ray fluxes. The morphol- 
ogy of this type of emission in such SNRs is expected to fol- 
low the structures of regions of enhanced density of target 
protons, not the structures in the SNRs where initial protons 
are accelerated. In SNRs which are bright in 7-rays, surface 
brightness distribution of proton-origin 7-ray emission may 
not therefore be expected to follow the radio and/or non- 
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thermal X-ray images of SNRs, like it is observed in IC443, 
W28 or CTB 37B. 

In contrast, the VHE 7-ray and hard X-ray morpholo- 
gies are observed to be well correlated in the cases of 
RX J1713.7-3946, Vela Jr., SN 1006 and possibly RCW86. 
Thus it could be that TeV 7-rays reflect the same struc- 
tures where the radio and nonthermal X-ray emission arise. 
In this scenario, electrons with energies of tens TeV may be 
responsible both for the (synchrotron) X-rays with energies 
of few keV and for the inverse- Compt on (IC) 7-rays with 
energies of few TeV which are observable by H.E.S.S. 

The H.E.S.S. image o f SN 1006 reported recently 
(Na umann-Godo et al.ll2009h reveals a very good correlation 
between X- and 7-ray maps. Can the correlation between IC 
7-ray and synchrotron X-ray images really be an argument 
for leptonic origin of VHE 7-ray emission? In the present 
paper, we make use of the spatially resolved analysis of the 
radio and X-ray data of SN 1006 to generate images with 
the possible appearance that this SNR would acquire if the 
whole TeV 7-emission were due to leptonic IC process. 

Since the purpose of the present analysis is to be as 
much model independent as possible, our work is mostly 
based on experimental results, without involving models of 
SNR dynamics, electron kinetics and evolution, etc., con- 
trary to what has been carried out in previous approaches 
to the problem (lRevnoldjll998l:[Fulbright fe Revnoldslll990l : 
iRevnoldsl 12004 lOrlando et all bo07l : IPetruk et al.1 l2Q09bF ). 
Radio and nonthermal X-ray emission contain information 
about the accelerated electrons, their distribution inside 
SNR, maximum energies, etc. The method we propose ex- 
tracts most of the important properties, which are needed 
to synthesize an IC 7-ray image, from the radio and X-ray 
data. The major exception is the magnetic field (MF). In the 
absence of observational information about it, we consider 
three cases of possible MF configurations. 



2 METHODOLOGY 

Let us assume that the energy spectrum of electrons holds 
the following relation 



N(E) = KE~ S exp(-£/£„ 



(1) 



where N(E) is the number of electrons per unit volume with 
arbitrary directions of motion, E the electron energy, K the 
normalization of the electron distribution, s the power law 
index and E max the maximum energy of electrons acceler- 
ated by the shock. This equation neglects small concave-up 
curvature of the spectrum predicted by efficient shock accel- 
eration but allows us to be in the fra mework of the met hod- 
ology of the X-ray spectral analysis ([Miceli et al 1 120091 . sr- 
cut model was used). The concavity results in a small bump 
around E max which leads mainly to some increase of the IC 
flux, which we are not interested in. It is not expected to af- 
fect the pattern of the gamma-ray brightness obtained with 
our method. Simpler spectrum N(E) = KE~ S is valid for 
the radio emission. The emissivity due to synchrotron or IC 
emission is 
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where p is the radiation power of a single electron with en- 
ergy E, e is the photon energy. The strength of magnetic 
field B is involved only in the synchrotron emission process. 

The simplest way to reach our goal is to use the delta- 
function approximation of the single-electron emissivities 
applied to spectrum Eq. (pQ). Namely, the special function 
F appeared in the theory of synchrotron radiation is substi- 
tuted with 
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where v is the frequency, v c (B,E) — c\BE 2 is the char- 
acteristic frequency, c± — 6.26 x 10 18 cgs. This results in 
synchrotron radio and X-ray emissivities with the following 
dependencies: 
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where z^break is 

Z^break = Vc(B , E n 
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With Eq. J4} and Eq. J5}, we approximate the relation be- 
tween the radio and X-ray synchrotron emissivities: 
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However, exponentially cut off electron distribution Eq. JTJ) 
convolved with the function approximation for the single- 
particle emissivity, Eq. J5}, underestimates the synchrotron 
flux from the same electron distribution convolved with the 
full single-particle em issivity, Eq. Q , at frequencies v > 
30z^break (see Fig. 3 in iRevnoldsl (|l998h , z^break is marked as 
zy m there). 

In this paper, we use zy x = 2.4 keV = 5.8 x 10 17 Hz. The 
range of v hre ^ in S N 1006 is found to be (0.06-=- 1) x 10 17 Hz 
(jMiceli et al. 2009). Thus, we are working with z/ x w (6 
100)z^break- Using Eq. (0 , we may therefore underestimate 
the real X-ray flux in ~ 10 times in regions where ^break 
i s sm all. Thus, as it is pointed out bv lRevnolds &; Keohand 
(1999), the approximate Eq. (0 is not robust at highest 
frequencies. 

We suggest therefore an empirical approximation of the 
numerically integrated synchrotron emissivity Eq. J2}, i.e. 
emissivity of the exponentially cut off electron distribution 
convolved with the full single-particle emissivity: 



-(s-l)/2 



exp 



(— ) 

V ^break / 



(8) 



where /3 X = 1.46 + 0.15(2 — s). This approximation is quite 
accurate. Its errors are less than 18% for s = 1.8 2.5 and 

Z/ x < 10 3 Z^break- 

We assume that VHE 7-ray emission from SN 1006 is 
due to IC process in the black-body photon field of the cos- 
mic microwave background. The convolution of the electron 
distribution Eq. JT} with the ^-fu nction approx imation for 
the single-particle IC emissivity ([Petrukl l2009h is also in- 
accurate to describe the 7-ray radiation of electrons with 
energies around E max . Therefore, like in the case of the syn- 
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Figure 1. Radio image of SN 1006 at A ~ 20 cm. The color scale 
is in units Jy beam -1 . The white contour denotes the boundary 
of SN 1006 (see text). 



chrotron radiation, we use an approximate formula for the 
IC emissivity, too. 

Let us consider IC emission at 1 TeV. We developed an 
approximation for the numerically integrated IC emissivity 
J2]) of electrons with the energy spectrum JT]) at 7-ray energy 
1 TeV: 
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where Ac = 15 for 2 < s < 2.5 and Ac = 15 + 2(2 - s) for 
1.8 < s < 2. The error of this approximation is less than 
25% for ^max > 0.3 TeV. This approximation accounts for 
the Klein-Nishina decline where necessary. 

The maximum energy is related to z^break by Eq. JB} : 
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where C\ — c x 1//2 . Substitution of Eq. J9} with this E n 
and K from Eq. J2} results in 
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This expression relates the radio emissivity and the IC 7-ray 
emissivity at 1 TeV with only one unknown, B. It may be 
used in the same cases where the srcut model is applicable, 
that is, if the spectrum of electrons may be approximated 
by Eq. JT]) with s assumed constant from the radio to X-ray 
emitting electrons. 

The idea of our method is represented by Eq. ([TT]) . 
Namely, this expression may be used in a small region 
("pixel") of a SNR projection in order to relate the surface 
brightness in radio band and IC 7-rays. Having this relation 
applied to all "pixels" , we may predict the main features of 
the 7-ray morphology of SNR originated in an IC process. 
This procedure 'converts' the radio image to an IC one. 

Another possibility is to start from the hard X-ray map 
and 'translate' it into the 7-ray image in a similar fashion. 
Namely, substitution Eq. ([TT]) with q r from Eq. JS} yields 



Figure 2. X-ray image of SN 1006 in 2 - 4.5 keV. The pixel size 
is 8". The color scale is in units MOS1 counts s _1 pixel -1 . 
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However, an X-ray image can be used only if it is dominated 
everywhere by the nonthermal emission. 

Eqs. and (|12[) relate emissivities of the uniform 

plasma. We use these equations to deal with surface bright- 
nesses that are superpositions of the local emissivities along 
the line of sight. Strictly speaking, this may be done only 
for the thin rim around SNR edge where plasma is approx- 
imately uniform along the line of sight. However, this ap- 
proach may also be extended to deeper regions of SNR pro- 
jection (see Appendix). Since X-ray limbs are (and 7-ray 
ones are expected to be) quite thin and close to the edge, 
our method is able to correctly determine the location of 
the bright limbs in the IC 7-ray image of SN 1006. In the 
interior of SNR projection, we consider B and z^break as 'ef- 
fective' values for a given 'pixel'. 

It is interesting to note that Eq. ([TT)) [or Eq. ([T2))] 
may be solved for the value of B. In this way, the method 
proposed here may be used for deriving the effective (line- 
of-sight averaged) MF pattern in SNR from its radio [or 
synchrotron X-ray] and IC 7-ray maps. The distribution of 
^break may be obtained from the radio (q T ) and synchrotron 
X-ray (q x ) images by solving Eq. JS}, without the need of 
the spatially resolved X-ray analysis. 



3 EXPERIMENTAL DATA AND MODELS OF 
MAGNETIC FIELD 

In order to make use of Eq. or Eq. (|12[) , one needs: i) an 
initial image, i.e. the map of distribution of the synchrotron 
radio (or X-ray) surface brightness, ii) the distribution of 
^break obtained from the spatially resolved spectral analy- 
sis of X-ray data and iii) the distribution of the effective 
magnetic field over the initial image. 
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Figure 3. Break frequency z^break : azimuthal profile (left) and image (right, the color scale is in units of Hz). 




Figure 4. Map of effective MF used in calculations for the model 
MF1. £ = 0.9. This map is rotated on 90° for the other scenario, 
MF2. The color scale is in units of B Q . 



We use the high resolution radio image of SN 1006 at 
A ~ 20 cm (Fig. [1]) , produced on the basis of archival Very 
Large Array data combi ned with Park es single-dish data pre- 
sented and described in llPetruk et al.l (|20Q9ah . 

As for the X-ray image, we use t he 2.0-4.5 keV XMM- 
Newton EPIC mosaic obtained by iMiceli et all (|2QQ9h . 
shown in Fig. [2] The mosaic has been produced by com- 
bining the 7 public archive XMM-Newton observations of 
SN1006, and using the PN, MOS1 and MOS2 cameras. The 
total effective EPIC exposure times ranges between 70 and 
140 MOSl-equivalent ks. 

A contour plotted on the images delineates the bound- 
ary of SN 1006. It corresponds to the level of 10% of the 
maximum brightness in the soft (0.5 — 0.8 keV) X-ray map 
of SNR (see Fig. 1 in lMiceli et al]l2QQ9h . 

The azimuthal profile (Fig. [3Jl) and image (Fig. [3J)) of 
the break frequency was derived on the basis of spatially re- 
sol ved X-ray spectral fitting results. The reader is referred 
to I Miceli et al.l (|20Q9h for the details on the procedure for 



X-ray data analysis and creation of this image. The same 
spectral fits also show that a, the radio spectral index, is 
between 0.47 and 0.53 everywhere around the shock. There- 
fore, we take s = 2a + 1 = 2 to be constant in SN 1006. 

We consider three models for magnetic field inside SNR. 
SN 1006 is rather symmetrical. In the procedure of MF map 
simulation, SNR is assumed to be spherical, with the radius 
equal to the average radius of SN 1006. For technical rea- 
sons, MF is fixed to the postshock value also outside the 
boundary of the spherical SNR (Fig. |4}. This allows us to 
deal correctly also with regions of SN 1006 which are larger 
than the average radius. 

Classical MHD description corresponds to the unmodi- 
fied shock theory. It takes into account the post- shock evo- 
lution of MF and the compression factor which increases 
with the shock obliquity. In this case, two possible orienta- 
tions of ISMF are considered. Namely, NW-SE (Fig. g]) in 
model MF1 (equatorial, or barrel-like, model) and NE-SW 
in MF2 (polar caps model). ISMF is assumed to be constant 
around SN 1006 with the strength B Q = 10 fiG. This value 
is choosen to give, in models MF1 and MF2, the postshock 
magnetic field 20 40 \iG, a value which follows from esti - 
mations for downstream MF strength (e.g. IVolk et al.l l2008) 
and reported 7-ray flux ([Naumann-Godo et al.ll2009f ). Such 
ISMF looks to be unrealistic at the position of SN 1006 far 
above the Galactic plane. A possibility to provide tens of \iG 
upstream of the shock would be the magnetic field amplifi- 
caion as an effect of efficient cosmic ray acceleration, which 
is out of the scope of this study. It should be noted however 
that the absolute value of the upstream field plays no role for 
the purpose of our paper. The aspect angle betw een ISMF 
and the line of sight is ta ken to be 70° for MF1 ([Reynolds! 
Il996l : IPetruk et al.ll2009ah and, for simplicity, also for MF2 
(we shall see below that the actual value of the aspect angle 
is not crucial for the purpose of the present paper, because 
even different models of MF lead to quite similar 7-ray pat- 
tern) . 

The procedure of generation of the average MF maps for 
MF1 and MF2 models is as follows. First we calculated nu- 
merically MHD model of Sedov SNR (jSedovll 19591 : iRevnoldsl 
Il998l : lKorobeinikovlll99lh . This gives us three-dimensional 
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Figure 5. Radio image of SN 1006 at A ~ 20 cm, from Fig. 03 
smoothed with Gaussian with 2' sigma. 



distribution of MF inside SNR. An effective MF in a given 
'pixel' of SNR projection is taken as a straight average for 
this 3-D MF distribution along the line of sight, account- 
ing for the azimuthal orientation and an aspect angle of the 
ambient MF in respect to the observer as well as the fact 
that most of emission arise right after the shock. Really, in 
order to generate map of effective MF one should calculate 
the emissivity-weighted average. Such an approach requires 
however the knowledge of 3-D distribution of emitting elec- 
trons within SNR which is unknown until one makes the 
full modelling from basic theoretical principles. In contrast, 
the scope of the present paper is to 'extract' the structure 
of radiating material from the observational data. There- 
fore, our procedure consists in approximate calculation of 
emissivity-weighted MF, without considering the 3-D distri- 
bution of relativistic electrons. In fact, most emission comes 
from a rather thin shell with thickness ~ 10% of SNR ra- 
dius. Therefore, in calculations of the average magnetic field, 
we consider only this part of SNR interior, namely the in- 
tegration along the line of sight is within regions from 
to i?, with £ < 1 (Fig. [4| . A choice of £ is rather arbitrary. 
It is apparent from calculations that contrasts between the 
outer regions and the interior of the IC image depend on 
this choice. The preference to the value of ( does not alter, 
however, the main features of the predicted IC morphology 
of SNR. Note that the azimuthal variation of the brightness 
is not affected at all by £ for radii of SNR projection > (R. 
In other words, the position of the limbs on the synthesized 
IC 7-ray image may not be altered by the particular choice 
of C- Nevertheless, in order to determine the most appropri- 
ate value of C, we made the full MHD simulations of Sedov 
SNR with model of evolu tion of the rela tivistic electrons 
in the SNR interior from iRevnoldsl (|l998h . Then the map 
of the emissivity-weighted average MF was produced from 
these simulations and compared with our maps of effective 
MF derived for different (. In this way, we found that the 
value C = 0.9 provides the good correspondence between MF 
maps in these two approaches. 

The third model of MF is relevant to the nonlinear ac- 
celeration theory with the time-dependent MF a mplification 
and the high level of turbulence (Bohm limit; IVolk et ah! 
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Figure 6. X-ray image of SN 1006 in 2-4.5 keV, from Fig. [2] 
smoothed with Gaussian with 2' sigma. 



2004). The quasi-parallel theory assumes in this case that 
the turbulence is produced ahead of the shock, not down- 
stream. The compression of the (already turbulent) mag- 
netic field then d o es not depend on the o riginal obliquity 
(IVolk et all 120031 : IBerezhko etHI l2QQ2h . iRakowski et~atl 
(2008) argue that shocks of different initial obliquity sub- 
ject to magnetic field amplification become perpendicular 
immediately upstream. ISMF is therefore assumed to in- 
crease on the shock by a large factor (due to compression 
and amplification), the same for any obliquity. In addition, 
in model MF3, we assum e MF to be approximate ly uniform 
everywhere insid e SNR (IBerezhko &; Volkl [2QQ4T ) . with the 
strength 150 fiG (jKsenofontov et al.ll2005f L 

Theoretical work on magnetic-field amplification start- 
ing with iBell &; Lucekl ([200 lh focuses on shocks which are 
originally parallel far upstream, with some implications 
that the process is less effective for (initially) perpendic- 
ular shocks. In this scenario, obliquity dependence of the 
post-shock MF would be opposite to the classical one. Re- 
ally, MF amplification is expected to follow acceleration ef- 
ficiency which decrease with the obliquity. If so, limbs in 
SN 1006 should correspond to the largest post-shock MF. 
Such MF morphology is qualitatively represented by MF1 
model (Fig. 13]). 

Some other notes to our methodology are in order. All 
initial maps were homogenized to the same size, orientation, 
resolution and pixel size. Eq. (|11|) is applied to each pixel 
of the initial images. The maximum brightness on images is 
fixed at the maximum value in the histogram distribution 
which has at least 10 pixels. The minimum is fixed at the 
level 1/100 of the maximum value. This is true for all the 
images, including the synthesized images, the observed X- 
ray and radio ones, and excluding the MF and z^break maps. 

Some images are smoothed to fit the resolution 
of H.E.S.S. The resol ution used is FWHM = 4.75' 
(Aharonian et al ] 12005(f ). so the Gaussian sigma is 2 . The 
role of smoothing is visible on Figs. [5] [6j to be compared 
with radio and X-ray images on Figs. [1] [2] 
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Figure 7. X-ray image of SN 1006 at 2.4 keV generated from the 
radio one (Fig. [T} with the use of Eq. JS) and smoothed to 2' 
Gaussian sigma. The color scale is normalized to the maximum 
brightness. 



4 RESULTS AND DISCUSSIONS 

To check our approach, we made use of Eq. JS} to generate 
the X-ray image of SN 1006 starting from the radio one. 
We stress that we do not need any assumption about MF 
configuration for this test, everything needed for radio to X- 
ray conversion may be taken from observations. Really, we 
need just the input radio image (Fig. [I]), and the distribu- 
tion of the break frequency (Fig. [3}. The resulting image is 
presented on Fig. [7] It shows good correlation with the hard 
X-ray observations (Fig. [6}, confirming that the proposed 
method works well. It restores the main properties of the 
observed hard X-ray image: essential decrease of the thick- 
ness of the two bright synchrotron limbs in X-rays comparing 
to the radio band (Fig. [5]), correct position of the limbs and 
negligible emission from the interior. Therefore, the method 
is reliable to be used for simulation of the 7-ray images of 

snr£|. 

Synthesized TeV 7-ray images of SN 1006 due to IC pro- 
cess are presented on Fig. [8] (model MF1, b arrel-like SNR 
in cla ssical MHD or polar caps in non-linear Bell Sz Lucekl 
(|200lh approach), Fig.[9](MF2, polar-caps in classical MHD) 
and Fiff. [TOl fMF3, u niform MF in the SNR interior for 
iBerezhko et all (|2002h non- linear model). Images presented 
on the left panels were obtained from the radio map as ini- 
tial model, while those shown in the right panels have the 
hard X-ray map as the starting point. Middle panels repre- 



1 The agreement between the observational and the synthesized 
X-ray maps is not perfect because the 'effective' resolution of the 
synthesized X-ray image cannot be better than the 'resolution' 
of the image of ^break- 

The effective resolution of ^break image 
is defined by the size of 30 rim regions used for spectral analysis 
(see Fig. 1 in lMiceli e t al. 2009) that in turn is determined by the 
photon statistics. Namely, the radial resolution in the image of 
^break is limited, since the radial profile of the cut-off frequency 
is assumed to be constant inside the rim regions. Note that the 
azimuthal 'resolution' is better than radial (see profile of ^break 
on Fig. 3). 



sent 'radio-origin' 7-images of the left panels smoothed to 
the resolution of H.E.S.S. 

Let us first consider the 7-ray morphologies obtained 
from the radio image. Two arcs dominate in all three MF 
configurations. Their locations correspond to limbs in radio 
and X-ray images. This confirms that correlation between 
TeV 7-ray and X-ray /radio morphologies may be consid- 
ered as direct evidence that the 7-ray emission of SN 1006 
observed by H.E.S.S. is leptonic in origin. 

Geometry of MF essentially different from those we 
considered might result in different predicted 7-ray images 
of SNR. Nevertheless, our results demonstrate that if MF 
strength varies within factor ~ 3 around the shock (Fig. 2} , 
any configuration of MF results in double-limb IC 7-ray im- 
age of SN 1006. 

For areas with the same radio surface brightness, higher 
MF implies lower IC 7-ray brightness, Eq. (|11[) . We have 
therefore different azimuthal contrasts (i.e. the ratio of max- 
imum to minimum brightness around the rim) in 7-ray im- 
ages for different models of MF. In particular, in model MF1 
the strength of the postshock MF is maximum in NE and 
SW regions (Fig. [4} where the radio brightness has maxima. 
At the same time, MF is smaller in faint NW and SE regions. 
This leads to small azimuthal contrast of brightness between 
bright and faint regions in 7-ray image (Fig. [8) left). In the 
opposite model MF2 the strength of the postshock MF is 
maximum in the faint NW and SE regions, that results in 
the largest brightness contrast (Fig. [9] left). Two arcs are 
therefore more pronounced in the 7-ray image for models 
MF2 and MF3. 

The differences in 7-images for three models of MF 
are not so prominent after smoothing them to the resolu- 
tion of H.E.S.S. In all three cases (Fig. l8l fT0l center panels), 
there are two bright limbs in the same locations correspond- 
ing to limbs on the smoothed X-ray map shown in Fig. [6] 
The synthesized images can be also directly compared t o 
the H.E.S.S. map of SN 1006 (|Naumann-Godo et alj|2009h . 
Good correlation between the synthesized and the observed 
images allows us to prefer leptonic origin of TeV 7-ray 
emission of this SNR. The uncertainties introduced in our 
method by the lack of knowledge of the real MF inside the 
remnant do not alter this correlation, as it is clearly seen 
from our synthesized images. 

VHE 7-ray images obtained from the initial hard X-ray 
map (Fig. IHl llOl right panels) are quite similar to those from 
the initial radio map, except for the configuration MF1. This 
fact reinforce the goodness of our method. The reason of 
differences between the 'radio-origin' and the 'X-ray-origin' 
7-ray images in the MF1 case (Fig. [Slight) is the contribu- 
tion of the thermal X-ray emission to the hard X-ray image 
which was used, in the SE and NW regions of SN 1006. 
Namely, our fitting show that the fraction of thermal emis- 
sion i n the overall 2-4.5 keV flux in the SE region is about 
50% (jMiceli et alJl2009h . The prominent but localized NW 
bright sp ot is completely dominated by the thermal X-ray 
emission (|Vink et al.ll2003h . So, the 7-ray brightness in these 
regions is overestimated in our synthesized images. This ef- 
fect is not prominent for MF2 and MF3 configurations be- 
cause MF is large enough in SE region to visually decrease 
the brightness there. Thus, the X-ray map may be used as 
initial one in our method only if it is completely dominated 
everywhere by the nonthermal emission. 
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Figure 8. Predicted IC morphology of SN 1006 at photons with energy ITeV, for the model MF1. IC image generated from the radio 
map (left); the same smoothed to 2' Gaussian sigma to fit the H.E.S.S. resolution (centre); IC image generated from X-ray map and 
smoothed to 2' Gaussian sigma (right). The color scales are normalized to the maximum brightness. 
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Figure 9. The same as on Fig. [8] for the model MF2. 
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Figure 10. The same as on Fig. [5] for the model MF3. 
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Figure 11. Azimuthal profiles of the observed radio (dashed line) 
and X-ray (crosses) compared with the synthesized 7-ray (solid 
lines) surface brightness (as from the input radio image), from 
images smoothed with 2' sigma Gaussian. Background galaxy at 
the azimuth around 90° is excluded. 



In order to quantify differences in images due to differ- 
ent models of MF, the azimuthal profiles from our synthe- 
sized 7-ray images are compared on Fig. [TT]with the corre- 
sponding azimuthal profiles derived from the VLA+Parkes 
radio image (Fig. [5} and the XMM-Newton X-ray image of 
SN 1006 (Fig. [6]). The profiles have been derived from maps 
smoothed to the HESS resolution. A large (from 8 / to 20' 
from center) annulus centered on the remnant is divided 
in 40 sectors. The values of brightness plotted is obtained 
by integration inside the sectors. Observational profiles are 
corrected for the background which is determined inside the 
central 0' — 7' circle. 

Fig. 1111 shows that the radio azimuthal profile may not 
be used as a tracer for the azimuthal variation of the IC 
7-ray brightness because it reflects only the changes in MF 
strength and density of relativist ic electrons, while X-ray 
profile accounts also for the variation of the electron maxi- 
mum energy which is important for 7-ray emission as well. 

Fig. HH re veals a very good match between positions and 
shape of the limbs in the hard X-ray map seen by XMM- 
Newton and in all synthesized 7-ray images. If the future 
reports of the H.E.S.S. collaboration reveal that the profile 
of VHE 7-ray brightness in SN 1006 closely match the XMM- 
Newton X-ray profile, then Fig. [11] could suggests that the 
preferred model for SN 1006 could be M F2, i.e. the bright 
limb s in this SNR might be polar caps ([Rothenflug et aD 
2004) in classical MHD description^]- However, it seems that 
statistical errors affecting modern H.E.S.S. data may not 
allow us to distinguish between the three models of magnetic 
field. We expect however that direct comparison of the X-ray 
and 7-ray observations with more sofisticated simulations 
(e.g. involving nonuniform ISMF) allows one to find a correct 
three-dimensional geometry for SN 1006 and ISMF around 
it. 



2 Note, that the polar-caps model does not explain the azimuthal 
profil es of the radio brightness, if ISMF is uniform (|Petruk et al.l 
2009a|). 



5 CONCLUSIONS 

We propose a general method to synthesize an IC 7-ray im- 
age of SNRs from their observed radio map and the results 
of spatially resolved X-ray spectral analysis, and we apply 
it to the supernova remnant SN 1006. The method is based 
on the fact that the surface brightness distribution of the 
synchrotron radio and X-ray emission of SNRs contains in- 
formation about the distribution and properties of acceler- 
ated electrons which are responsible for the 7-ray emission 
as well. We have derived analytical expression to calculate 
the 7-ray image from radio or X-ray surface brightness map, 
Eq. (|lip . It is used for each pixel of the input image, knowing 
the surface distribution of z^break and assuming some config- 
uration of magnetic field. 

X-ray and TeV 7-ray photons are radiated by electrons 
with almost the same energies. It seems therefore to be more 
natural to use the hard X-ray map as input rather than the 
radio one. However, we show that the contribution from the 
thermal component in the initial X-ray image may result in 
incorrect prediction for the IC brightness distribution. Thus, 
the hard X-ray map may be used only if it is completely 
dominated everywhere by the nonthermal emission, which 
is not the case of the XMM-Newton 2.0-4.5 keV image of 
SN 1006. 

On the other hand, the usage of the radio map as input 
may introduce some errors due to possible change of the 
spectral index in the electron spectrum. The applicability 
of the method may be tested by generating a synchrotron 
X-ray image from the radio one with Eq. JS}. If synthesized 
distribution of X-ray surface brightness correlates well with 
the observed one (as in the case of SN 1006) then both the 
srcut model in XSPEC and our method for IC 7-ray image 
may be used for a given SNR. 

Our synthesized VHE IC 7-ray images of SN 1006 are 
quite similar to that r eported in the publication of the 
H.E.S.S. collaboration (jNaumann-Godo et all l2Q09h . This 
fact favours a leptonic scenario for the TeV 7-ray emission 
of this SNR. However a hadronic origin cannot be ruled out 
in view of the measure d ISM densities, consist ent with a 
hadronic scenario (e.g. iKsenofontov et aD l2005h . If this is 
the case, the observed TeV brightness map may reflect the 
distribution of protons with energies > 2.4 TeV which inter- 
act with compressed ISM downstream of the shock. 

The present spatial resolution achieved by H.E.S.S., 
prevents us to ultimately disentangle between the three con- 
sidered configurations of MF in SN 1006: in all three cases, 
the two arcs on our 7-ray images are in the same location 
as in radio and X-rays. If MF strength varies no more than 
factor ~ 3 or so (as on Fig. [4} around the shocks, IC 7- 
ray image of SN 1006 should have two limbs located in the 
same regions as in X-ray map, independently of the actual 
azimuthal configuration of MF. The reason, in accordance 
with Eq. (fTT]) , is the contrasts in (radio or X-ray) brightness 
and ^break which dominate any moderate azimuthal varia- 
tion of MF. The use of the observed ratios of the radio sur- 
face brightness and the break frequency between NE and SE 
regions in Eq. (|11[) shows that a variation of MF Bne / ^se 
larger than a factor of 4 may reverse the location of bright 
IC limbs. 
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APPENDIX A: APPROXIMATE DESCRIPTION 
OF SURFACE BRIGHTNESS 



Here, the approach developed in IPetruk et al. I (|2009al lbh is 
adopted for general description of the surface brightness. 
The surface brightness of a spherical SNR projection at dis- 
tance g from the center and at azimuth ip is 



S(Q,<p) 



f R 

\ I q(a, Bo 

Ja(g) 



rr a da 



(Al) 



where q is emissivity, o = Oo(^ 5 r / 4>o) is the shock obliq- 
uity, c/) an aspect angle, r and a are Eulerian and Lagrangian 
coordinates, r a the derivative of r(a) in respect to a. The 
emissivity in synchrotron or IC process is 



dEN(E)p(E,v) 



(A2) 



In the ^-function approximation of the single-electron emis- 
sivity p(E, v\ we may write that 



q oc N(E m )B* 



(A3) 



where E m is an energy of electron which gives maximum 
contribution to radiation at a given frequency zy, x = 1/2 
for synchrotron and x = for IC emission. 

Energy spectrum of electrons N(E) evolves in a differ- 
ent way downstream of the shocks with different obliquity, 
i.e. N(E m ) = N(E m ;a,S ). In Sedov SNR, this evolution 
may approximately be expressed by the two independent 
terms (for details see lPetruk et al.|[2009al lbh 



N(E m] a, Go) « N a (a)N@(g, B eff ) 



(A4) 



where O e ff = S ((p,l, 4> ). The similar relation holds for 
MF: B(a, O ) ~ B s (0 o )B a (a) where B s is the immediately 
post-shock value. This allows Eq. (|A1[) to be written as 

S(g,<p) ocN@(g,e eS )B s (e eS ) x 



I N a (a)B a (a)*- 

Ja(o) 



rr a da 



(A5) 



'a(g) \/r 2 - Q 2 

where integral depends on g only. In other notations, 

S(q, (p) « g e ff (q, (f) • T(q) (A6) 

where X is an integral in (|A5[) devided by N a (g). The accu- 
racy of this approximate formula increases toward the edge 
of SNR projection where the bright limbs we are interested 
in are located. 

It is important that the factor I does not depend on hjxa ; 
but only on the coordinate in the projection. This means 
that X is almost the same in a given position of the radio, 
X- and 7- ray images and we may use Eqs. (181). (pJ3 and (H2l) 
written for emissivities in order to relate surface brightnesses 
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in each 'pixel'. The factor X is also independent of ip in this 
approximation. Therefore, the azimuthal variations of the 
surface brightness S at a given g may be directly represented 
by the azimuthal variations of the effective emissivities. This 
provides justification for discussion in Sect. 2] However, the 
radial contrasts in brightness should account for the radial 
changes in T which is unknown until one considers detailed 
3-D MHD model of SNR. 



